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Bose-Einstein condensation (BEC) of composite diquarks in quark matter (the color supercon- 
ductor phase) is discussed using the quasi-chemical equilibrium theory at a relatively low density 
region near the deconfinement phase transition, where dynamical quark-pair fluctuations are as- 
sumed to be described as bosonic degrees of freedom (diquarks). A general formulation is given for 
the diquark formation and particle-antiparticle pair-creation processes in the relativistic framework, 
and some interesting properties are shown, which are characteristic for the relativistic many-body 
system. Behaviors of transition temperature and phase diagram of the quark-diquark matter are 
generally presented in model parameter space, and their asymptotic behaviors are also discussed. 
As an application to the color superconductivity, the transition temperatures and the quark and 
diquark density profiles are calculated in case with constituent /current quarks, where the diquark is 
in bound/resonant state. We obtained Tc ~ 60 — 80 MeV for constituent quarks and Tc ~ 130 MeV 
for current quarks at a moderate density (pt ~ 3po)- The method is also developed to include 
interdiquark interactions into the quasi-chemical equilibrium theory within a mean-field approxi- 
mation, and it is found that a possible repulsive diquark-diquark interaction lowers the transition 
temperature by ~ 50 %. 

PACS numbers: 12.38.-t, 12.38. Mh, 12.38.Lg, 03.75.Nt 
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I. INTRODUCTION 
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Several arguments based on the quantum chrome-dynamics (QCD) show that the hadronic matter makes 
transition from the hadron phase to the quark-gluon one at sufficiently high temperature and/or large 
baryon-number density. In particular, recent intensive studies on the quark matter suggest that, in low 
temperature, it becomes a condensed state of quark pairs called the color superconductor (CSC) , where 
the color SU(3) symmetry is spontaneously broken 0,0, In model calculations, the most attractive 
quark-quark interaction induced by gluon exchange appears in the total spin singlet (J = 0) and the 
color anti-triplet (3) channel, so that the CSC phase is considered to be a state where the quark pairs 
with these quantum numbers are condensed coherently in a similar mechanism with the BCS theory Q . 

Generally, there exist two limiting cases in the fcrmion-pair condensation: weak- and strong-coupling 
regions in the strength of the two-fermion correlation. In the weak-coupling region, which appears in 
the high-density quark matter because of the asymptotic- freedom of QCD j5j, the coherence length £ 
becomes larger than the interquark distance d, and the quark matter becomes the weakly-interacting 
fermion system where the mean-field picture are applicable. The strong-coupling region appears in 
the moderate- or low-density quark matter (especially in the neighbor of the deconfinement transition 
point). Large pair-fluctuations develop [|| and the coherence length is suggested to be very small; £/d ~ 
0(1 — 10), even around the moderate-density region [~ O(10)po] from the analysis of the improved ladder 
approximation 7J. This implies that, for low-density region [~ O(l)po], £/d -C 1 and the strong quark- 
pair correlations are described by bosonic degrees of freedom, which we call "diquark" in this paper. 
Thus, for the description of this diquark region of the quark matter, the incorporation of two-particle 
correlations beyond mean- field levels is essentially important. 

Several methods have been proposed to incorporate the correlation effects. A typical one is that of 
Nozieres and Schmitt-Rink, where thermodynamic potentials are calculated in the normal phase using 
the in-medium ladder approximation, and the fermion number density is shown to include a contribution 
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from the bosonic part through the two-body in- medium scatterings |8|, |9j . The transition temperature 
is obtained from the zero-frequency pole of the two-particle correlation function (Thouless criterion for 
second-order phase transitions). 

The quasi-chemical equilibrium theory is an another approach developed by Schafroth, Butler, and 
Blatt [l0( for the fermion-pair condensates. Using a cluster expansion technique for the partition function, 
the fermion-pair condensates are shown to be in chemical equilibrium states between the isolated fermions 
and the bosonic difermions identified with the in-medium resonant/bound states of two fermions. This 
method was firstly invented to describe the experimental phenomena of electric superconductivity of 
metals in 1940', before the birth of BCS theory In modern sense, this quasi-chemical equilibrium picture 
corresponds to the strong coupling superconductor < 1), and the transition temperature, which is 
obtained as the BEC-Tc of difermions under the equilibrium, is shown to be equivalent with that in 
Nozicrcs-Schmitt-Rink theory in the strong-coupling limit (£,/d <C 1). 

In the strong-coupling region of the quark matter, we can expect that the similar quark-pair fluctuations 
develop and the CSC is described by the diquark condensate in equilibrium with isolated quarks. 

The diquark correlation effects have been proposed by Jaffe and Low 0, in relation with the 
phenomenology of the low-energy hadron-hadron collisions and the multi-hadron spectroscop y. Also, 
roles of composite diquarks in the quark-gluon plasma have been studied by Ekelin, et al. |23l |24| . 
especially for the hadronization process in the high-energy e + e~ collision; they considered the space-time 
developments of the quark-gluon plasma as an adiabatic process, and discussed the possibility of the 
diquark BEC rather at high-T stages in that process. Also the diquark correlations have been discussed 
to study the stability of diquark matter but without regard to diquark BEC [25l I2H I2H l28l l29l | . 

In this paper, we make a general formulation of the quasi-chemical equilibrium theory for the relativistic 
many-fermion system, and apply it to the strong-coupling quark matter with composite diquarks. Our 
special interest is on diquark BEC that should be the CSC in the strong-coupling region, and we clarify 
the CSC states as diquark BEC at finite temperature and density within the quasi-cquilibrium theory 
in less model-dependent way. In Sec. II, we formulate the quasi-chemical equilibrium theory for the 
difcrmion formation and the particle-antiparticle pair-creation processes in the relativistic many-body 
system. We also show some interesting properties for the BEC of difermions which are characteristic in 
the relativistic system. In Sec. Ill, applying the method to the quark-diquark matter, we show Tq of 
the diquark BEC, a phase diagram on quark/diquark masses, and asymptotic behaviors of them. We 
also discuss Tq of the color superconductor using phenomenological values of parameters. In Sec. IV, 
effects of the interdiquark interactions are discussed in a mean-field approximation. Sec. V is devoted to 
summary and outlook. Throughout this paper, we use the natural units: c = h = ks = 1- 



II. RELATIVISTIC QUASI-CHEMICAL EQUILIBRIUM THEORY 

A. Molecular formation/annihilation process 

To develop the quasi-chemical equilibrium theory, we consider molecular formation/annihilation process 
in a fermionic matter: 

F 1 +F 2 ^(F 1 F 2 )=B, (1) 

where F\^ are fermions with masses mp = rriF 1 — itif 2 , an d B is a composite boson of F\ and F 2 
with mass mg; the composite boson B should be a bound state (mj < 2m p) or a resonance state 
(tub > 2mp). We call it "molecule" in both cases. 

Here we take the quasi-particle picture for the fermionic matter where the system consists of quasi- 
fermions (F\ and F 2 ) and quasi-bosons (B). The two-body interactions between fermions bring about 
the two-body correlations in the matter; their major effects are for the binding/resonance energy of the 
composite boson B and the residual interactions among quasi-particles. This quasi-particle interactions 
are generally regarded to be weak, and we neglect it in the first part of this paper. In Sec. IV, we 
introduce the quasi-particle interactions and discuss their effects within a mean-field approximation. 

The equilibrium condition for the process Q is given by 

M-Fi + VF 2 = MB, (2) 

where /j-f 1 ,f 2 ,b are the chemical potentials of the particles F±, F 2 and B, which should be functions of 
the temperature T and the particle-number densities np 1 ^p 2 _B- For free uniform gases, they are given by 



the Bose/Fermi statistics: 
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where /3 = 1/T. Using the relativistic energy-momentum dispersion relations: e a = W m 2 + k 2 (a 
Fi,F 2 ,B), eqs. Ilol 11) become 
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n F a = 7TT / /^-„„ -i , T de = m F F R (-Pn Fet ,m Fa P), (a =1,2) (6) 
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In the above equations, the relativistic Bose/Fermi functions and _Fr are defined by 
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where we use scaled variables /? = m(3 and e = e/m, and the Juttner variable 9 that is defined by 

e = coshe nsmj. 

The boson chemical potential /^s satisfies //^ < ms from the positive norm condition of ns, and 
the boson density |J5J has a singularity at this upper bound, where a phase transition occurs to the 
Bose-Einstein condensate. In the BEC region (fiB = tub), the boson density consists of the condensed 
and thermal parts, rig = + n^ h \ which are given by = m B Bp(—f3mB,mBf3) and n^P — 

riB-n^M- 

In the process Q, the particle-number conservations for Ft 2 give constraints for the densities: 

n Fl +n B = n Fld , n F2 + n B = np 2:t , (9) 

where n Fl 2 .t are the total number densities of -Fi,2i which consist of isolated fermions and those included 
in the composite boson B. 

Solving eqs. J5J and (0 with Ij5l6|l . we can obtain the densities riB,F lt F 2 in equilibrium at temperature 

T. 



B. Particle-antiparticle pair-creation/annihilation process 

The energy scale in the quark-diquark system is so large that the antiparticle degrees of freedom are 
generally important. Here we develop the quasi-chemical equilibrium theory for the particle-antiparticle 
system. 

The fermion-antifermion (quark-antiquark) or boson-antiboson (diquark-antidiquark) pairs can be cre- 
ated through the radiation process mediated by the gauge boson (gluon); 

G^P + P, {P = F 1 ,F 2 ,B) (10) 

Thus the equilibrium condition for the above process is given by fip + up = uq, where u P p G are the 
chemical potentials of the particle P, the antiparticle P and the gluon G. We put uq = because the 
gluon is a massless particle, so that the chemical potentials of P and P have a relation: up = —fip. 
In the process (|1(JH . the difference between the particle and antiparticle densities should be conserved: 



nb,p=np-n P , (P = F 1 ,F 2 ,B) 



(11) 
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where n&.p corresponds to the baryon-number density of quark or diquark degrees of freedom. As given 
in (l-'iKil) . the boson/fermion number-densities of the particle and antiparticle at the temperature T are 
given by 

n B = m B B R (-(3n B ,m B 0), n B = w? B B R {^\x B ,m B f3), (12) 
n F a = Tnp a F R (—p(j, Fa ,rn Fa P), n Fa = m Fa F R (f3[i Fa ,m Fa {3), (a = 1,2) (13) 

where B R and F R have been defined in J7JSJ). 

Eq. (JSJ) shows that the boson chemical potentials are found to have the upper bound; fi B B < m B . 
From the relation fi B = — [J-b-i they are also found to have the lower bound: 

- m B < H B ,B < m B, (14) 

where \i B > (fi B < 0) corresponds to rib, B = n B — n B > (nb,B < 0). Furthermore, the upper and 
lower bounds of /is in (|14|l correspond to the boson BEC (fi B = Wb) and antiboson BEC (/i^ = m B ) 
respectively, which means that bosons and antibosons cannot be condensed to the BEC state at the same 
time in equilibrium states. 

The transition temperature Tc of the boson BEC (the antiboson BEC) is determined by with 
[i-B = m B (fJ-B = ~m B )- 

nb,B = Kb - n B } fJiB=±mB ^ T=Tc . (15) 



C. Molecular formation/annihilation processes of many kinds of particles 

Here we consider a general system of many kinds of fermions F a and molecular bosons Bp which 
are bound/resonant states of two fermions. The fermions are assumed to be transmuted through some 
processes; F a «-* F a i . As is same with the previous simple case, the total number density of fermions 
n Fjt is to be conserved in molecular formation/annihilation process of Bp, and is given by summing all 
species of fermions n Fclt t, which also amounts to summation of isolated fermions n Fa and those included 
in composite bosons n B „: 

n F , t = n Fa>t = n Fa + 2}^ n Bg ■ (16) 

a a p 

To obtain the equilibrium condition for the above condition, we take the Hclmholtz free-energy density 
of the system: 

/ = F/V = f T ,v + ^2 Mf q n Fa + fJ'Bp n Bfj 

a j) 

-K \ n Ftt -Y,n Fa -2J2 UB I 1171 

a p 

where /t,v is the density-independent part, and the last term has been added to expose the constraint 
(|16|) with the Lagrange multiplier K. Differentiating eq. I|17f) with respect to n Fa and n B[jl one obtains 
for V V ^ 

^M-Fq = 2if = fi B/3 . (18) 

This means that the chemical potentials of fermions and bosons take the same values /i F and \x B in an 
equilibrium state, which obey the condition J2Jl; 2fi F — [i B . 

The range of the chemical potential is fi B = fx Bg < m Bfj for any Bp and, without loss of generality, 
the boson B\ is assumed to be the one with the smallest mass (most-stable composite); m Bl < m B/3 ^ 1 . 
Then, the range of \i B becomes 

Hb{= 2^f) < m Bl . (19) 

It shows that, at the upper bound of the chemical potential ((J, Bl = fJ. B = m Bl ), the BEC occurs for the 
boson B\, but there exists no BEC singularity for the other bosons because Mb^i = Ms = m B 1 <m Bf3 . 
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In summary, the BEC occurs only for the boson with the smallest mass below the transition temperature 
Tc, and the other bosons with larger masses (less-stable composites) cannot be condensed to the BEC 
states at any temperature and density (one-BEC theorem). 

The critical temperature Tc for the Si's BEC can be determined from the condition at T = Tc: 



n Ft . 



R 

Below Tc, the density of B\ consists of the thermal and condensed parts: 

A th ) 3 



L B 1 



m Bl B R (-(3m Bl ,m Bl P), 



n F ,t 
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(20) 

(21) 
(22) 



where np,t — Y) a n F a ,t is the total fermion number, and n Bfs ^F a are given by I|12I13|) again. 

In the case of the system with antiparticle degrees of freedom, antibosons cannot be condensed to BEC 
state if the baryon-number density is positive. This can be also understood by one-BEC theorem. The 
ground state contributions of bosons {B{) and antibosons {B\) are explicitly given by 
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(23) 
(24) 



which means that the antiboson is effectively heavier than boson by 2fi B (> 0) and BEC occurs only for 
bosons Bi by using one-BEC theorem (In case that the baryon-number density is negative, y, B < 0, only 
Bi can then be condensed to BEC state at the lower bound jj. B = —m Bl ). 



III. DIQUARK BOSE-EINSTEIN CONDENSATION IN QCD 



In this section, we apply the relativistic quasi-chemical equilibrium theory developed in the previous 
section to the color superconductivity of quark matter in BCS-BEC crossover regime [ll|,[l2l> where quark- 
pair fluctuations develop dynamically to create diquarks (composite bos ons) . A possibility of quark-quark 
correlations has also been suggested in low-energy hadron physics [l3l Il4j. and its extrapolation to the 
deconfinement phase can be considered to become diquark degrees of freedom in dense region [2^, l28l . . 



A. Applicability of quasi-chemical equilibrium theory for quark matter 

Before going to actual calculations we should evaluate the applicability of quasi-chemical equilibrium 
theory for interacting fermionic systems, which was originally discussed by Schafroth, et al. with respect 
to electron gas in superconducting metals |lf| : 

1) particles of compositions (fermions F and composite molecules B) should be quasi-particles, 

2) correlations between FB and BB are unimportant, 

3) higher-order correlations are not stronger than the pair-correlations considered, 

4) any localizations to liquid droplets with large volume exclusion (gas-liquid transition) do not occur 
before the formation of quasi-molecules, 

5) composites sufficiently obey the Bose statistics as expected with small size of composites relative 
to the interparticle distance. 
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Now we discuss the condition 1) for quark matter. The lifetime of a particle relates to various decay 
processes in medium, and to recognize quark as quasi-particle its lifetime should be at least longer than 
inverse of the QCD energy scale, that is, 0(1) fm. The lifetime of quarks tends to be shorter near the 
Fermi surface, which can be roughly estimated from the imaginary part of quark self-energy, e.g., |l5j : 
t ~ 30 _1 MeV -1 . In our model these quarks with short lifetime near Fermi surface due to strong pair 
fluctuations are described as bosonic molecules, and residual quarks in the deeper interior of the Fermi 
sphere have explicitly longer lifetime than the QCD scale. Thus, quarks (antiquarks) being outside of 
pairs due to medium effects can be identified as quasi-particles. As for the diquark composites, its lifetime 
as quasi-particle is expected to be longer than that of quarks in crossover regime, as Shuryak and Zahed 
suggested the existence of various binary bound states in the whole regions just above the deconfinement 
temperature because QCD coupling strength is found not to be so small as one has believed so far [16). 
Thus, diquarks should also be identified in medium. 

Concerning other conditions 2) and 3), higher-order correlations rather than pair ones are not treated in 
this Sec. IIL This treatment corresponds to the Gaussian-type approximations like Nozieres and Schmitt- 
Rink approach 0, 0, 0] . The effects of higher-order correlations, especially interdiquark interactions, 
will be discussed in Sec. IV with sincere generalization of quasi-chemical equilibrium theory. Condition 
4) is also discussed in the end of Sec. IV. 

With regard to condition 5), the density dependence of coherence length £, characterising typical size 
of a Cooper pair, relative to interquark distance d is estimated by using the improved ladder approxima- 
tion 0, which suggests that £/d < 1 for pi, < 10po- This means that (/del for lower density region 
Pb ~ 3po concerned, and our quasi-chemical equilibrium approach with bosonic modes is expected to 
work well for these strong couplig regimes of quark matter. 

The critical coupling strength for the formation of diquark bound states is also calculated by using 
effective four- Fermi interaction model with Gaussian approximation [j^. and critical value G/(n 2 /A 2 ) ~ 
0.92 (A is ultraviolet cutoff scale) and its relavance for real quark matter is suggested therein. Because 
of smooth crossover, this critical value of coupling strength corresponding to the unitary limit in cold 
atomic systems brings no sudden changes for physical quantities. The essential quantity for the behavior 
of BEC-BCS crossover is always the ratio of the coherence length to the interparticle distance. 

Although above quantitative argumets come from some specific model calculations, recent analyses on 
QCD experiments, e.g., relativistic heavy ion collision 01 an d lattice QCD 0], also suggest that even in 
deconfinement phase QCD holds strong coupling at almost whole region of the phase diagram relevant to 
possible observations. We expect that our present approach works well in such strong coupling regimes 
in more general way, especially around the region close to the phase transition line where correlations 
between quarks should develop largely and the strong coupling picture becomes essential for the pairing 
problem. In this context three-quark correlations in addition to pairing will be important issue around 
this region, which can also be treated in our approach. 

Here we comment on the BEC-BCS crossover in beta-equilibriated charge-neutral quark matter, where 
mismatches of Fermi surfaces exist in different flavors. Even in the weak coupling regime, various phases 
have been proposed in that matter in relation with the chromomagnetic instability problem: the gapless 
superconductor |l9|. the crystalline color superconductor [2(J, the gluon condensed state [5l|, and so 
on. The strong coupling limit has also been discussed of the two-flavor quark matter with Fermi-surface 
mismatch [2^ ■ These considerations suggest that the Cooper-pair formation should still be favored in the 
case if the energy gain by the pairing exceeds the energy differences of the mismatched Fermi surfaces. 
Thus we can expect the crossover into the diquark BEC in the quark matter with finite mismatches 
for sufficiently large coupling strength and, at the temperature around Tc, the mismatch effects are 
expected to be small with large pairing energy. In low temperature region, however, the Fermi-surface 
mismatches makes the unpaired fermions remain with some flavor, and they should form new Fermi 
surfaces; the calculation of the profile changes of them around the crossover region should be a very 
interesting problem in future. 

In the present paper, we take the two-flavor quark matter with same Fermi surfaces to discuss typical 
features of diquark BEC; however, our approach can be extended to the case of the finite Fermi- surface 
mismatches and their effects in the diquark BEC and crossover will be discussed in another paper. 

B. Baryon-number conservation and chemical equilibrium condition 

In what follows, we take u,d-quark matter around the moderate density region in QCD, and consider 
the chemical equilibrium between quarks and composite diquarks in the 2SC-type pairing state (color 
antisymmetric 3, flavor antisymmetric, and total-spin J = 0) Q Q The other types of diquarks 
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can also be created in the quark matter using less attractive channels, e.g., diquarks in color-3, flavor 
symmetric and the J = 1 pairing state [33l l34| . But, according to the one-BEC theorem derived in the 
previous section, the diquark BEC can occur only for the smallest mass (most stable) diquark, i.e. in the 
2SC-diquark state of u and d quarks, so that we consider only the 2SC-diquarks. Generally the diquarks 
which do not participate in the BEC would have the effect of lowering the critical temperature, through 
the total-number conservation, in comparison with the case of one kind of diquark. We also assume that 
the quark and diquark masses do not depend on the temperature and density, but the present formalism 
is also applicable to the case with effective masses dependent on these parameters. 

In the quark-diquark matter, the baryon- number density pb(= rib/3) corresponds to the conserved total 
fcrmion- number density nj? )f in Ijllfl or lltj|) . where rib is the regularized baryon- number density in quark 
degrees. In the chemical equilibrium between quarks and 2SC-diquarks, it is given by 

n b = 3 • 2 • 2 • {n q - n q ) + 2 • 3 • 1 • 1 • (n d - n 3 ) , (25) 

where n q . q are the number densities of quark and antiquark, and n d d are those of diquark and antidiquark. 
The numerical factors in 1251) come from the degeneracy of the internal degrees of freedom: the color- 
triplet (3), the flavor-doublet (2), the spin-doublet (2) for the quark-density term, and the compositencss 
(2), the color-anti-triplet (3), the flavor-singlet (1), the spin-singlet (1) for the diquark-density term. 
Eqs. H12I13|) gives the densities n qjq and n d d in equilibrium at temperature T: 

n q = m q F R (-/3fi q ,m q l3), n q = m q F R ((3u q ,m q l3), (26) 
n d = m d B R (-l3fi d ,m d f3), n d = m d B R (Pp d ,m d /3), (27) 

where m q ^ d are a quark/ antiquark and a diquark/antidiquark masses, and pp = —pp (P = q, d) has 
been used for the chemical potentials of antiquark and antidiquark. When the diquark (antidiquark) 
BEC occurs for T < Tc, the diquark (antidiquark) density consists of the condensed and thermal parts: 
nd = n d +n d (n d = +ni*^), which are obtained in the same way as eqs. H21I22|> . Especially these 

condensates and n§' are the superfluid density of CSC. If the baryon-number density is positive, 
diquark BEC can only be expected as mentioned in the previous section. 

The chemical potentials p q and p d should obey the quark-diquark equilibrium condition: 

2fx q = (28) 

which corresponds to J2J. The set of eqs. (|25I28|) determines a composition of the baryon-number density 
in the quark-diquark quasi-chemical equilibrium system. 



C. Transition temperature of quark-diquark system 

As shown in the previous section, diquark BEC occurs when the diquark chemical potential reaches 
the upper limit u d — m d . In the quark-diquark system considered here, the conditions for Tq are 

2^ 9 = Pd = m d , (29) 
n b = h ■ 2 • 2-(n q - n q ) + 2 • 3 • 1 • ~n ( , th) )} , (30) 

a j t=T c 

where we have used eqs. (|25I28|) . For any fixed values of the quark and diquark masses m q . d , and for the 
baryon-number density rib in (|30|l . Tc is obtained by solving these equations. 

Fig. d shows dependence of Tq on 2m q and m d for the baryon-number densities p b = 2po,3po,Apo 
where po ~ 1.3 x 10 6 MeV 3 is the normal nuclear density. These figures show that Tc increases with 
the baryon-number density; especially it can be read from the shrinkage of no-BEC region (area in the 
front side in the (2m 9 ,m £ i)-plane). This effect can be understood that Tc goes up with increase of the 
diquark density together with the baryon density, as seen in the behavior of Tc in the ideal bose gases 
(see Appendix lAl. 

These diagrams also show that the large diquark mass has an effect of lowering Tc- It is consistent 
with the results for one-component boson-antiboson gas (see Appendix B); Tc oc m^ 1 in nonrelativistic 

limit i|A8|l and cx m^^ 2 in ultrarelativistic limit i|A17(l . The increase of Tc for large quark-masses is 
attributed to the large diquark density brought by the increase in the binding energy (m d — 2m q ). 

In the very small mass region of diquark (m d ~ 0), which corresponds to the extremely tight binding 
case (strong-coupling limit), we find rapid increases of Tc, which is divergent at the massless limit m d — * 



8 




FIG. 1: Transition temperatures Tc of diquark BEC on (2m q , m,j)-plane of quark/diquark masses for the baryon- 
number densities pb = 2po,3po,4po (po is the normal nuclear density). 



(Fig.[2Ji). It can be understood from the range (|14|l and the condition (|28|l of the chemical potentials, from 
which we obtain pL q = [id — in the present limit, meaning that quark and thermal diquark components 



of the baryon- number densities vanish (n q — n q = 0, 



<th) 



= 0) and the conserved baryon-numbcr 



density is fully compensated by the condensed diquarks (n& = 2 • 3 • n^) at any temperature (that is, 
Tc — > oo). This peculiar feature is characteristic of the massless boson with the conserved quantity which 
is not satisfied by photons, et al.. Quark mass dependences of Tq are shown in Fig. |2p, and we can find 
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FIG. 2: Transition temperature Tc of diquark BEC: diquark-mass (mj) dependence (a) and quark-mass (m q ) 
dependence (b), at pb = 3po- Several lines are for different quark masses m q = ~ 750 MeV (a), and for diquark 
masses m d = 400 ~ 1000 MeV (b). 



that, for large values of m qi Tc converges into a finite limit that depends on diquark mass md'- Tc of free 
bose gas with mass m^. It is because the quark degrees of freedom arc suppressed in the limit of large 
quark mass. 
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D. Mass phase diagram of quark-diquark system 

In the previous subsection, numerical results (FigQ show that there exists the region in (2m q ,rrid)- 
planc where no BEC occurs. In Fig.0 we plot Tc — lines in [2m q , md)-plane for several baryon- number 
densities Pb(= rib/ 3). This line is the border of the BEC and no-BEC phases (below and above the line), 
and we call it mass phase diagram, where BEC phase is defined with finite Tc- 

no-BEC phase 

1000 
800 
> 600 

§ 

| 400 
200 


200 400 600 800 1000 
2m q (MeV) 

FIG. 3: Mass phase diagram of diquark BEC for pt, = (1 ~ 5)po- Curved lines correspond to critical borders 
between BEC and no-BEC phases. Dashed line is for 2m q — rrid, separating region 1 (bound state case ; 
rrid < 2m q ) and region 2 (resonant state case ; rrid > 2m,). 




Now we show an analytical derivation of the border in the mass phase diagram by taking the Tc 
limit in the formula. At T — 0, all antiparticles vanish, so that eq. (|25|l becomes 



n b = [3 • 2 • 2 • n q + 2 • 3 • 1 • 1 • rid] 



T^0 ' 



(31) 



n d th \ Just at 



The diquark density generally consists of the condensed and thermal part rid = 
T = Tc, the chemical potential reaches its upper limit pd — rrid but there is still no BEC component, and 
thermal diquarks are completely frozen in the limit of T — 0, so that we obtain rid = at T = Tc = 0. 
The quark density is a monotonously-increasing function of p qi and it takes the maximum value at 
T = T c = 0; < lax with fi q = p d /2 = m d /2. Thus eq. ffl} becomes 



From eq. ®), 



n b = 3 • 2 • 2 • rf 



the n™ ax 



is evaluated to be 



(32) 



1 

2^2 

1 

6^ 



:^/e 2 — m q ■ 6(p q — e)de 



H q =m d /2 



2m, 



3/2 



x 9(rrid — 2m q ). 



(33) 



In the case of m d < 2m q (region 1 in Fig.|3J), the above equation gives n" lax = (no quarks), which leads 
to the contradictory result of rib — with l|32() . It means that Tc > and all quarks are combined into 
diquarks with manifest advantage of binding energy and should be condensed to BEC state at T = 
(the BEC phase). 
When rrid > 2m q (region 2 in Fig. QJ, eqs. (|32I33|I gives 



nb = — j 



3/2 



(34) 



For the fixed value of rib, it gives the hyperbolic curve: the border between the BEC and no-BEC phases on 
(2m q , m^-plane as shown in Fig. [3] The behavior of this border is decided by the balance between kinetic 
energy loss of isolated quarks (Pauli-blocking effect) and resonance energy loss of composite diquarks. 
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E. Diquark condensations for constituent and current quarks: transition temperature 

Now we take a look at proposed values of the quark and diquark masses in quark matter, given by 
Anselmino, et al. |35| . where two kinds of the quark/diquark masses have been presented: 



case 1 : m d = 500 MeV, m q = 200 - 300 MeV, 



case 2 : m d = 200 MeV, 



0- 10 MeV. 



(35) 



The case 1 corresponds to dynamical chiral symmetry breaking which remains even in the deconfinement 
phase, and thus the large quark mass comes from the, as it were, constituent quarks. Whereas the case 2 
is for the chiral restored quark matter with small current quark mass. For the diquark masses, it should 
be noticed that the diquark is bound (m d < 2m q ) or resonant (m d > 2m g ) states depending on the value 



of 



in case 1, but is only the resonant state in case 2. In this paper, we take these two cases for the 



quark and diquark masses. 

Fig. 01 shows dependence of Tc on the baryon-number density in the cases 1 and 2. The larger quark 
masses get Tc higher because the diquark binding energy (ma — 2m ? ) becomes larger. This tendency is 
less apparent in Fig.^D compared with Fig.^, which is understood that the scale of current quark mass 
is too small relative to the energy scale of diquark BEC. 



(a) case 1. ( m rf =500MeV ) 



(b) case 2. ( m,,=200MeV ) 
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FIG. 4: Baryon-number density (pb) dependences of Tc for case 1 (a) and case 2 (b). 



In Fig. [5] comparison for Tc in the nonrelativistic and ultrarelativistic limits is shown in the cases 1 
and 2. The analytical representations of Tc in these limits are given by 



(a) case 1. ( m rf =500MeV, m,,=300MeV ) 



(b) case 2. ( m rf =200MeV, m,,= 10MeV ) 
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FIG. 5: Baryon-number density (pt.) dependences of Tc for case 1 with (m,d,m q ) — (500, 300) MeV (a) and case 
2 with (md,m q ) = (200, 10) MeV (b), accompanying with Tc in nonrelativistic and ultrarelativistic limits. 



n 2/3 

, 3.313-^, (nonrelativistic limit) , , 

T c ~ { m A— V ^ (36) 

1.732 /^_. (ultrarelativistic limit) 



The derivation of these formulas are given in Appendix A. 
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In Fig. [3^ (case 1), Tc is found to be approximated by the nonrelativistic formula in lower density 
region < 2po), and approaches to the ultrarelativistic one in the high density region; in the inter- 
mediate region, the gradual shift is shown between these two limits. It is a typical behavior of the 
relativistic system. On the other hand, in case 2, the small current quark mass makes the ultrarelativistic 
approximation valid from the very small density region (Fig. [5Jd). 



F. Diquark condensations for constituent and current quarks: density profiles of quarks and 

diquarks 

Here we discuss the temperature dependences of particle and antiparticle densities shown in Fig. 
(case 1) and Fig.[7|(case 2) at /?{, = 3po- The baryon-number density of quark is defined by p^ A = p q — p q , 

and that of diquark consists of the thermal part p£j = p d th ' — p^ h ' and the condensed part p^\ = p^j/ 
(Here we consider only the case of pf, > 0, so that no condensed antidiquarks exist; p^ — 0). 



(a) case 1. ( m £ ,=500MeV, m,,=300MeV ) (b) case 1. ( m rf =500MeV, m,=300MeV ) 




_q 5 I .... i .... i .... i .... i .... i .... i .... i .... I _i I .... i .... i .... i .... i .... i .... i .... i ... . 
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FIG. 6: T-dependence of the baryon-number densities (a) and the particle/antiparticle- number densities (b) for 
constituent quarks (case 1) at pt = 3po- Lines q and q are for p q ,q, and those D and D are for p 1* t ■ Lines [g — q] 
and [D — D] are for pb, q and ■ Line is for the condensed part p£\ = p^\ 



(a) case 2. ( m,,=200MeV, m,= 10MeV ) 



(b) case 2. ( m (/ =200MeV, m ( ,= 10MeV ) 
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FIG. 7: T-dependence of the baryon-number densities (a) and the particle/antiparticle- number densities (b) for 
current quarks (case 2) at p b = 3po- For conventions, see the caption of Fig. El 



In Figs.|Hland[7| p d begins to be finite at a sufficiently low temperature showing BEC phase transition. 
It should be noticed in particular that the case 2 is diquark BEC in the unstable resonance state of the 
current quarks (|35fl . which shows the stabilization of the resonance states in the matter through the 
equilibrium processes. The existence of the BEC is also confirmed in the saturation behavior of the 
quark chemical potential (Fig. |SJ; it shows that p q {= Md/2) reaches the upper values rrid/2 at Tc and 
keeps this value in T < Te- 
la case 1, the diquark is in the bound state (md < 2m g ) and the association into the diquark (q + q — * d) 
gains the binding energy, which gives ph, q = at T — > 0; all quarks are combined into diquarks. In case 
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FIG. 8: Temperature dependences of quark chemical potentials in case 1 and 2. 



2, the diquark formation increases the total energy of the system because the diquark is a resonance 
state (rrid > 2m q ), and increase of the quark Fermi-energy (and also the entropy effect) competes with it. 
Balancing these effects, the system does not so much favor the pairing into bosonic molecules, and quarks 
still remain even at T ~ (Fig, Ek)- These behaviors around zero temperature is generally discussed by 
using mass phase diagram (Fig. E| in the previous subsection. 

Now we turn to discussion for the high-T region. Taking case 1 first (Fig. Eh), we find that the thermal 
diquark density p d decreases from T ~ 100 MeV ~ (2m q — mj), around which the dissociation should 
be effective, up to T ~ 200 MeV, but at this temperature the particle-antiparticle pair creations begin 
to occur, and then the diquark (and any other particle) density increases again. In case 2 (Fig. 0d), we 
also find the density increase that can be attributed to the pair creation effects. These diquark-density 
increases are very characteristic phenomena in the relativistic condensates. Really as shown in Fig. El in 
case of no antiparticle degrees of freedom, all composite molecules should dissociate into isolated fermions 
due to entropy enhancement. 



(a) case 1. ( m rf =500MeV, m,,=300MeV ) 



(b) case 1 . 
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FIG. 9: Temperature dependences of the particle-number density (a) and quark chemical potential (b) in case of 
no antiparticle degrees of freedom. (md,m q ) — (500, 300) MeV, pb = 3po- For conventions, see Fig. El 



Next we analyze asymptotic behaviors of the quark and diquark densities at the high-temperature limit. 
In Fig. 110b and b, the high-T behaviors are shown of the quark/diquark baryon-number densities in the 
case 1 and 2. From these figures, we can find a feature common in both cases: the ratio pb,q/ Pb.d — ► 1/2 
at T — > oo. From the analysis given in Appendix B, we obtain a universal representation at the high- 
temperature limit: 



R 



q/d 



= f\q 
Pb.d 



(37) 



9c ■ Pq 



9st 




= 1A 



(38) 
(39) 
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(a) case 1. ( m £ ,=500MeV, m (/ =300MeV ) 



(b) case 2. ( m (/ =200MeV, m (; =10MeV ) 
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FIG. 10: Asymptotic behaviors of the quark/diquark baryon-number densities at large T and pb = 3po for case 1 
(a) and 2 (b). For conventions, see Fig. [5] 



where g q .d are the multiplicity factors for quark and diquark in the present system, which represent the 
degenerate degrees of freedom (color x flavor x spin). g c is the composite factor indicating that a diquark is 
constructed by g c x quarks, and g st is the statistical factor which reflects the statistical difference between 
quarks (fermions) and diquarks (bosons). In the present case, g q = 3 • 2 • 2, <jy = 3 • 1 • 1, g c = 2, and 
g st — 1/2. With these values the ratio l|37|l becomes R q /d — 1/2, which shows that the numerically 
obtained limit is universal, being independent of masses m q> d- 

As shown in Appendix B, the quantum statistics and the antiparticle degrees of freedom are essential in 
the derivation of the correct asymptotic limit l|37[l . Especially existence of the statistical factor g st in (|37() 
implies that the quantum statistics still remains even in the high-temperature limit. This is understood 
that particles and antiparticles are drastically produced through the pair creation with large thermal 
energy, and Boltzmann region corresponding to high-temperature and dilute limits almost disappears. 

Finally, we should mention consistency of the present method with real QCD at high temperature and 
density. Diquarks are dynamically-produced composite states by interquark interactions, which should be 
suppressed due to the asymptotic freedom in the high temperature and/or density. Furthermore, the Pauli 
blocking effect on the quark-quark interaction prevents quark pairs from forming composite diquarks, and 
the diquarks should become correlated pairs with a coherence length comparable with the interdiquark 
distance at the high density. Such effects of increa sing density for pairing mechanism is referred as 'Mott 
effect' discussed with nuclear matter in refs. 1^3, l38j . Since nuclear force has a finite repulsive core in 
general, nuclear matter should undergo superfluid-normal phase transition at a specific density with short 
interparticle distance comparable to the range of the hard core, after the BEC-BCS crossover as density 
increases. QCD does not have such a specific scale of interaction length in the deconfinement phase, 
thus one can expect that the system is still in the superfluid phase even with infinite density at least 
around zero temperature. In this high density regime, however, diquarks as composite degrees become 
unessential with large coherence length as discussed above. Thus applicability of the present method 
should be restricted in relatively low-density region, and incorporation of these effects in some fashion is 
needed for extension to higher densities. 



IV. EFFECTS OF INTERDIQUARK INTERACTIONS FOR DIQUARK BEC 

The quarks and diquarks discussed in this paper should be quasiparticles in quark matter; the inter- 
action effects are partially included in the existence itself of the composite states (diquark) and in the 
effective masses of the quasiparticles (quark and diquark). These effects have been included in the calcu- 
lations in the previous subsections, so that they are not for the free quarks in the real meaning. However, 
we have still a missing effect, which are also not included in the Gaussian- type approximations |3,l2(' 
that the quasiparticles can scatter into the different states through the residual interactions. In this 
subsection, we consider the diquark-diquark interaction within a mean-field approximation, and figure 
out its effects for the transition temperature of the diquark BEC. 

As shown later, in the single-component boson systems, the boson-boson interaction effect does not 
change Tc from that of the free bose gas l|A13|) within the mean-field approximation I3! j; i t is consistent 
with the observed Tc of atomic BECs, which is rather a strongly interacting system |4(|. The boson- 
boson correlation effects somewhat shift the value of Tc in the observed A-transition of the liquid He 4 : 



14 



T\ = 2.17K (eq. (jA"T3|) gives T c = 3.1 K), which is also shown in theoretical estimations |3al4ll4al43.l44|. 
In the diquark BEC, differently from the single-component boson system, the interaction effects are found 
to be effective for Tq still in the mean-field approximation through the chemical equilibrium process 
(q + q <-> d), so that the calculations in this level give a first-order approximation. 

We start with a model of interacting color-3 scalar fields <f> a (a = 1, 2, 3) [4lTl42T|. which represents the 
diquark degree of freedom: 



(40) 



where (/> 4 -term gives a diquark-diquark scattering effect (the Gross-Pitaevski approach). Here we have 
assumed that the higher-order scattering terms ((f) ' '"•) are renormalized into the two-body interac- 
tion 1451. It should be noted that the Ginzburg-Landau free-energy for the CSC phase also has a similar 
form |4* a l47 | . 

In the mean-field approximation, the interaction energy of the system is given by the expectation value 
of the interaction term in l|4U|) : 



{Hi) ee (A / d 3 x4(x)0 Q (x)0t( x )^( x) ) 





d 3 x(^(x)0 Q 


A 
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w 


E 

k,l,a,/3 





+ 



2ngnf 



'k"i 



2n£ 



2nl 



(41) 



where and n k are the occupation numbers of the diquark and the antidiquark states with a color 
index a and a wave vector k, and u>^_ = -\/k 2 + m 2 d . 

Now we consider the cases 1 and 2 in (|35[) just above Tc- In both cases, the antiquark densities can 
be neglected because Tq (a few tens MeV as in Figs. 0J is smaller than the diquark pair-creation energy 
(2md = 400 MeV ~ 1 GeV) . We also keep the leading-order term for the particle- number density because 
we consider the high-density region pB| . Thus the interaction energy (|41|) becomes 



(H T 



A 

W 

A 
4 



k,l,a,/3 W k W l 



V 



d 3 \td 3 \ n k n\ 
(2ir) 6 LOk^>\ 



+ 



d 3 \i 2n nk 
(27r) 3 m d LUk 



1 

V 



no 



(42) 



where rik = ^2 a and in the last line we have taken the continuum limit. It should be noted that, to 
manage the BEC singularity, we picked up the k = term no from the k integrals in the continuum limit 
in Eventually, the total energy of the system is given by 



E d = V 



d 3 k 

(2tt) 3 



cj k (n k + Uk) + (Hi), 



Differentiating with respect to we obtain the single-particle energy: 



D 



(2^) 3 SE d 
V 5n(k) 

d 3 k 



D 



«k 



^k 2 + m 2 d 
n 



(27r)3 



m d V 



(43) 

(44) 
(45) 



To evaluate (|45p. we concentrate on case of a small coupling constant and a large diquark mass; B < m 2 d . 
Generally we can say that this condition should be satisfied for small value of A, but its exact evaluation 
is not easy so that we take it as a working hypothesis in the present paper. Under the condition B < m 2 , 
eq. (|44() has its minimum at k = 0. More general arguments including wide range of A and m d will be 
made elsewhere. 

As discussed in the previous subsection, for the sufficiently larger diquark mass m d , the nonrelativistic 
approximation should be valid at the temperature close to Tc- In such a case, the occupation number 
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rik (the Bose distribution function) should vanish for k 3> TTI4, so that we can put ^/k 2 + m d ~ nrid in 
the first term of (|45|l in good approximation. Thus, we obtain the result: 



£k — ► Wk 2 + m 



i 2 + 

d + 2m 2 



d 3 k 

\3 



"0 



(2tt) 3 * y 

k 2 +^ + ^|3(nf» (46) 

where n d th ^ and are the thermal and condensed densities of diquark. It should be noted that the 
antidiquark density nk has been dropped in (|46() because it is less effective than that of diquark in the 
nonrelativistic approximation. 

In the mean- field approximation, the single-particle energy appears, e.g., in the distribution functions 
as £k — Mtii where the chemical potential fid is introduced as a Lagrange multiplier to the conserved 
baryon- number Q = J d 3 xjo(x) with the baryon- number current j M = i{4>*d^(j) ~ 4>d^4>*) f° r the diquark 
degree of freedom. Using eq. ljl6*|l. we obtain 



v^^+ ^3 („fJ+ 4 0) ) -/* 



= ^k 2 +m 2 -^, (47) 

where the first term corresponds to the free-particle dispersion and the effective chemical potential \x* d is 
defined by 



A 



2m d 



*3 (»<*>+ »»?>), (48) 



which is a procedure of the chemical-potential renormalization j39j . 
With eq. i|47[l . the boson distribution function 127f) becomes 

n d = m 3 d B R (-f3fx d , m d [3). (49) 

In this case, Tq is determined by the condition fi d = md in (|49|) : = m d Bn(—f3md, rrid/3), so that the 
formula for Tc itself with fixed nd is completely the same as in the non-interacting bose gas; therefore 
Tq of single bose gas should not be changed at all through the chemical-potential renormalization [39j 

Since we are considering the multi-component system with the chemical equilibrium (q + q d), the 
equilibrium condition 128fl must hold for the bare chemical potentials fi q ,d- 

2 N = ^ = ^ + ^3 U d th) + 4 0) ) • (50) 
zm d v ' 

In the equilibrium, this equation and the density constraint ll'ol) should be solved with the Bose and 
Fermi distribution, l|49f) and l|26[). The transition temperature Tq is obtained by putting further condition 
li*d = m d- It is clear that, differently from the single-component case, the interaction term in l|5l)[l gives 
a A-dependence in Tq, which is the unique property of equilibrium system. 

If the interaction is repulsive (A > 0), it has an effect of decreasing the diquark density because the 
interaction term (|46J) gives the energy increase in the large diquark density, so that Tc gets lower with 
coupling strength A. This is clearly found in Fig. 1111 where the numerical results are given for the 
A-dependence of Tc with the repulsive interaction in the cases 1 and 2 in l|35|l . 

The diquark-diquark interaction has been discussed for the color-3 diquarks by Donoghue and 
Sateesh |25|, where the value of A = 60 ~ 100 is estimated from the mass difference of the nucleon 
and A with the assumption of the quark-hadron continuity |25j . The positiveness of A is also suggested 
in the P-matrix method 0] . 

Figs. IT2*1 shows phase diagrams for the interacting diquark system with the suggested values A = 60 ~ 
100. From these results, we can expect that residual diquark-diquark interactions have the effect to lower 
Tc by ~ 50 % in comparison with that in the non-interacting cases. 

Here we should comment that the other interactions except the diquark-diquark one might give im- 
portant effects, e.g., a possible quark-diquark interaction should generate baryonic collective modes in 
relation with three-quark correlations. As for these higher order correlations, we can discuss these effects 
by extending our method, referring to higher-order cluster expansion technique [To| . 
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FIG. 11: Interdiquark-interaction dependence of Tc in case 1 (a) and 2 (b) for the several baryon-number densities. 
Coupling constant A is defined in 1401 . 




Finally we should discuss a possibility of the transition of quark matter into classical liquid with large 
volume exclusion. Schafroth, et al. 10] pointed out that the quasi-chemical equilibrium theory cannot 
be applied to the interacting fermionic system which undergoes the condensation into the classical liquid 
before the formation of difermion molecules. To evaluate it for quark matter we now estimate the de 
Boer parameter A* [4^ which gives the ratio of the quantum energy fluctuation to the potential depth e: 

A'= * =-^, (51) 
ay/me aVmc 2 e 

where a is a range of interaction and m is a particle mass. Large A* corresponds to large quantum 
fluctuation relative to the potential depth e. Taking typical values of several parameters for quark 
matter; a ~ 0(l)fm, mc 2 ~ O(10)MeV, e ~ O(10)MeV, A* is roughly estimated to be O(10 2 ) which 
is very large relative to that of the quantum liquid He (A* ~ 2.64). This large quantum fluctuation 
implies that quark matter cannot be condensed into classical liquid and it must be a quantum gas or, 
rather, quantum liquid. As qualitative feature of liquid He 4 can be well described by free bose gas model 
calculations, we expect, in this sense, the quasi-chemical equilibrium picture should be applicable for 
quark matter. 



V. SUMMARY AND OUTLOOK 



We have generally developed the quasi-chemical equilibrium theory for the system of isolated fermions 
and composite bosons of fermion pairs in the relativistic framework, and derived the conditions for 
occurrence of the composite-boson BEC. Especially a general property of BEC in a multi-component 
system is suggested; only the lightest composite bosons can be condensed to the BEC state (one-BEC 
theorem). 
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An application has been done for the diquark condensates in quark matter at moderate densities, 
and only 2SC-type diquarks are treated by virtue of the one-BEC theorem. In this framework we have 
discussed the color superconductivity as the diquark BEC in a possible BCS-BEC crossover regime. 
We found the mass phase diagrams which determine the occurrence of BEC at a certain temperature 
for various mass values of quarks and diquarks. With phenomcnological values of model parameters, 
we obtain Tq< as a few tens MeV, which is lower than those by model calculations within mean-field 
approximations, ~ O(100)MeV. In any case, Tq has a large dependence on quark and diquark masses, 
and it varies from to oo along with — > (strong coupling limit). 

At the high-T limit, antiparticles produced through the pair-creation processes suppress the dissociation 
of diquarks into quark pairs, and the ratio of the number densities are found to approach a universal 
value R q /d — 1/2, which reflects that quantum statistics still remains even in the high-T limit. 

We have also discussed effects of the diquark-diquark interaction, and found within the mean-field 
approximation that a possible repulsive interaction lowers Tc by ~ 50 % in comparison with that of 
non-interacting case through a chemical-potential renormalization. 

The other correlation effects beyond the mean-field approximation, e.g., quark-diquark and three-quark 
correlations, should also be important in the BCS-BEC crossover problem. Combining the method of the 
Beth-Uhlenbeck approach with the present quasi-chemical equilibrium theory, we should discuss the 
correlation effects and the crossover problem from less model-dependent point of view. The study along 
this line is now ongoing and will be presented near future. 
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In this appendix, we derive some formulas on the critical temperature of the BEC transition for 
single bose gas in the three limiting cases: 1) nonrelativistic limit, 2) ultrarelativistic limit with a mass 
correction, and 3) ultrarelativistic limit with negative-charge particles (antibosons). 

Here we consider the function Br defined in (JJJ , the denominator of which can be expanded as 
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APPENDIX A: ANALYSIS OF BEC-TRANSITION TEMPERATURE 




(Al) 



where [ib = Ms /m. The integral in the above expansion is evaluated to be 




(A2) 



where K v {z) is the modified Bessel function: 




(A3) 



and we have used a formula: K v ^\(z) — K u +i{z) 
We obtain 



{2v/z)K v (z) HH. Substituting eq. (H) into (JUJ, 




(A4) 
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The transition temperature (3c = tub jTc is obtained from eq. (JSJ) with p,B — /tub = 1, so that 



ml ^ e i0c 



n B = m%B R (-p c , W = r4V -~-K 2 (0 c ). (A5) 

Z7T 1=1 L Pc 



Using the asymptotic expansion of K v (z) for z — + oo 49]: 

-k 



K * z) V2/ £r („-* + §) ft! ' (A6) 



we obtain the asymptotic expansion of (| A5|) for large (3c (small Tc): 

i « r(f + *)c(fc^ 3 



where At = y2Tr/3/m B is the thermal de Broglie wave length and the ^-function is defined by £(z) = 
y^ ^ n l~ z ■ Taking the first term in the right-hand side of l|A7(l . we obtain the nonrelativistic limit for 
TcM- 

n B ^Tc ~ C(3/2) ~ 2.612. (A8) 

To obtain the ultrarelativistic limit, we use the asymptotic expansion of K n (z) with positive integer n 
for z — > 0: 



1 , (tn — h — \\\ , z \ 2fc-n 

^(^Ec- 1 ) - fc , - Q) • ( A9 ) 



Taking the first term of the above expansion and using it for l|A5[) . we obtain the density formula for 
small (3c ■ 



n *~^Y,7WT 3 = "T^Mc), (A10) 

n l^i WcY k 2 P c 



where B a {v) is the Appell function 

"1 = 

T(a)J e*+ 



1 f°° z a - x dz ^ e- lu 

B " M = rwl :^-t=£— < A11 » 



Using the expansion B 3 (—(3c) ~ C(3) + C(2)/?c for /3c — > (high- Tc limit) in eq. ijAlOfl . we obtain the 
ultrarelativistic limit for Tc'. 



n B 



03) 02) 



= ^T c + m B ^Tg, (A12) 

where C(2) = tt 2 /6 ~ 1.6449 and ((3) ~ 1.2020. It should be noted that the first term in the above 
equation does not include tub', in the ultrarelativistic limit, the boson mass is neglected and the energy- 
momentum dispersion relation is approximated to be e ~ |k|. The second term in eq. (|A12|) corresponds 
to the small mass correction. 

In summary, we obtain the formulas of Tc of the relativistic BEC: 

3.313-^-, (nonrelativistic limit) (A13) 
2.017n]/ 3 . (ultrarelativistic limit) 
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Finally we consider the boson-antiboson system, where the similar calculation with l|15|) gives the 
expansion for baryon-number density : 

m| ^ sinh(l/3/2 B ) ~ 
n b , B = n B - n s = — > -~ K 2 {lj3). A14 

The formula for Tq is obtained by putting (1b = 1 and (3 = (3c- Using eq. (|A6fl . we obtain the asymptotic 
expansion for large (3c'- 



nb ' B = ^^ P^)h{ C { k + l) - B *+*/2Wc)}(2Pc)- k , (A15) 



where -Bfc+3/2, defined by (|A11J) . corresponds to the antiboson contributions. In the nonrelativistic limit 
(low Tc), Bk +3 / 2 (2(3c) can be neglected because B a (z) ~ e~ 2 , so that we obtain the same formula 
as eqs. (|A7IA8ll where ns is replaced with n b .B- This means that the antiboson degree of freedom is 
suppressed and we can put n bt B ~ ns in this limit. 

In the ultrarelativistic limit (small (3c) , using l|A9|) . eq. (|A14I) becomes 



2m r w „ z/r 2 „ 
^C(2)|| = ^— - 



- ^C(2)^ = ^C(2)/iBl", (A16) 



where sinhz/z ~ 1 has been used. Substituting \±b — 1 and /3 = (3c, we obtain the ultrarelativistic 
formula for Tc'. 



2m|C(2) 9 C(2), 
2 c n 



n b , B ~ _a" ^, = 2rn B ^YT^. (A17) 



Comparison with eq. I|A12|) shows that the leading terms from n B g cancels out so that the result depend 
on tub in spite of the ultrarelativistic limit. 

In summary, we obtain the formulas of Tc for the boson-antiboson system: 

^2/3 

, 3.313—^, (nonrelativistic limit) 
T c ~ { m A— (A18) 

1.732 (ultrarelativistic limit) 



APPENDIX B: ASYMPTOTIC FORM OF DENSITY PROFILES 



First we calculate the asymptotic expansion of fermion-number density; calculations go in the same 
way as for the boson density formula (|A14|I : 



, sinh(l/3fip) 



n b , F =n F -np = ^ V(-l)'- 1 y -^^K 2 (l(3), (Bl) 

^ ^ 1(3 



' f , ,/- 1 

1=1 

where the extra factor (— 1)' _1 comes from the factor +1 in the denominator of the fermion-density 
JSJ reflecting Fermi statistics (Bose statistics gives the extra factor (+1)' _1 = +1 which is embedded 
in eq. QA14|l l. Expanding the function K 2 and taking the leading-order term in (|A9(I . we obtain the 
fermion-density formula for small (3 : 

2mK(2)MF _ 2 C(2) 2 

n »,F — ^2 ~ 2~ ' (B2) 

where IX^ -1 )'" 1 /^ = i 1 - 2 1_a )C(a) has been used to give the factor C(2)/2 in (|B2)l . In c ase of 
Bose statistics, factor /l a = C(«) from t he analy sis of eq. (|AT4|l gives C(2) in JXTBJ. The 

difference between the factors £(2) and C(2)/2 in eqs. I|A16IB2| tells that the effect of quantum statistics 
still remains even in high-T limit, which is discussed in eq. (|37|l . 



20 



Now we apply the boson and fermion density- formulas for high-T region, l|A16|) and (|B2|) . to quark- 
diquark system, and evaluate the high-T limit of the baryon-number density ratio. Using eqs. (|A16IB2|I . 
the baryon-density ratio of isolated quarks to diquarks is given by 



R, 



q/d 



Pb,q 
Pb,d 



(l/3)g q n b>q 



{2/3)g d n btd 



9 st, 



(B3) 



9c ■ (J>q 
9st 



E 



(-1) 



l-l 



I 2 



(+1) 



i-i 



i=i 



P 



1/2, 



(B4) 
(B5) 



where 1/3, 2/3 are the baryon numbers of quark and diquark, and g q ^ are the multiplicity factors for quark 
and diquark in the present system, which represent the degenerate degrees freedom (color x flavor x spin). 
g c is the composite factor indicating that a diquark is constructed by g c x quarks (in the present case, 
g c = 2), and g s t is the statistical factor which reflects the statistical difference between quarks (fermions) 
and diquarks (bosons). 
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